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Abstract 

Multiplicity correlations between the current and target regions of the 
Breit frame in deep-inelastic scattering processes are studied. It is shown 
that the correlations are sensitive to the first-order perturbative QCD effects 
and can be used to extract the behavior of the boson-gluon fusion rates as a 
function of the Bjorken variable. The behavior of the correlations is derived 
analytically and analyzed using a Monte Carlo simulation. 
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1 Introduction 



Long-range correlations in rapidity have been studied for many years in e + e~, /i + p, 
up, up, pp, pp, 7r ± p, K + p, and e + p collisions in forward-backward event hemispheres. 
Results for e + e~, /i + p processes indicated that such correlations are small at energies 
studied in At LEP1 energies, DELPHI, OPAL and ALEPH collaborations 

observed positive long-range correlations, mainly due to heavy quark pair production 
For up and up processes, the correlations are rather small and negative ||. 
For pp, pp [|7| and 7r ± p, K + p || collisions, the correlations are positive and increase 
with y/s. Recently, it was shown that the long-range correlations defined in the 7*P 
center-of-mass system of diffractive e + p collisions are positive ||. 

In this paper we discuss the measurements of the long-range correlations in 
neutral current deep-inelastic scattering (DIS) e + p collisions using the Breit frame 
jnj. We estimate analytically the correlations from the first-order QCD (Sect. |3|) 
and compare them with a Monte Carlo simulation (Sect. f|). 

2 Definitions 

The event kinematics of the deep-inelastic processes is determined by the 4-momentum 
transfer Q 2 = —q 2 , and the Bjorken scaling variable x = Q 2 /(2Pq), where P is the 
4-momentum of the proton. To study the correlations, we use the Breit frame. In 
the quark-parton model (QPM), the Breit frame provides the maximum separation 
between the radiation from the outgoing struck quark and the proton remnant. In 
this frame the incident quark carries momentum Q/2 in the positive z-direction 
and the outgoing struck quark carries Q/2 in the negative ^-direction . The phase 
space of the event can be divided into two regions. All particles with negative p z 
components of momenta form the current region, which is analogous to a single 
hemisphere of e + e~ collisions. In the QPM, all these particles are produced due to 
the hadronization of the struck quark. Particles with positive p z are assigned to the 
target region, which is associated with the proton remnant (see Fig. 1). 

Long-range correlations deal with the problem of a possible interdependence 
of different, well-separated phase-space regions of multiparticle production. For 
the Breit frame, it is natural to ask whether the current and target regions are 
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independent of each other. Assuming that a color flow between the struck quark and 
the proton remnant at the fragmentation stage cannot produce strong correlations, 
one may think that the radiation of the struck quark should be independent of the 
target region. Below we shall show that, generally, such an expectation is not correct 
when QCD corrections to the QPM are considered. 

To study a possible correlation between the current and target region, we use 
the covariance: 

p = cov(n c , n t ) = (n c n t ) - (n c ) (n t ) (1) 

where n c (nt) is the number of particles in the current (target) region, (. . .) is average 
over all events. 

There are a few well established statistical properties of the p: 

1. Fully independent particle production in the current and target regions implies 
p = 0. 

2. Since p represents a degree of linear stochastic dependence between n c and 
n t , more complex forms of interdependence are not described by this variable. 
Therefore, the fact that p = is still not evidence for absence of correlations. 

3. Positive correlation leads to p > 0, negative one produces p < 0. 

The forward-backward correlations have been studied earlier in terms of the 
parameter b which is given by the slope in the linear relation (rip) = a + bn^, where 

(^b) is the number of particles in the forward (backward) hemisphere. The 
parameter b is directly related to the covariance p and can be defined |T]J using the 



standard deviation op of the multiplicity distribution in the forward hemisphere, 

b = of p. (2) 

For the present study we shall use p rather than b. While the statistical content 
of these quantities is same, the covariance is easier to obtain both analytically and 
experimentally. 

Having established all these notations and definitions, below we shall analytically 
estimate the covariance p using first-order QCD effects. 
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3 QCD expectation 



Going from the quark-parton model to the first-order QCD description means to con- 
sider more complex processes. In Born approximation, Boson-Gluon Fusion (BGF) 
and QCD-Compton scattering (QCDC) can significantly contribute to the overall 
topology of the events (see Fig. 2). Both processes lead to (2+1) partons, where 
"1" denotes the proton remnant. While in the QPM the quark has p z = Q/2 before 
and p z = —Q/2 after the collision, the first-order QCD processes make the collision 
in the Breit frame no longer collinear, although the target and current regions are 
still well defined operationally. 

Following ||12|| , different event configurations of (2+1) jets in the Breit frame 
are shown in Fig. 3. The configurations are defined according to the longitudinal 
momenta of jets induced by the first-order QCD processes. For the first configuration 
both jets from the first-order QCD processes are moving to the current region. For 
the next two topologies (see b) and c)), the two jets are produced back-to-back with 
one jet in the current and one in the target region. (Both jets may have different 
longitudinal momenta, but this does not affect the QCD estimate of correlations to 
be made below.) For the topology d), the two jets are produced in the target region. 
Note that, for the Breit frame, the latter topology is not in conflict with longitudinal 
momentum conservation [|1^ . 



Let now estimate the covariance p analytically. If h particles migrate to the 
target region, one can rewrite ([!]) as 

p = ((n c - h) (n t + h)) - (n c - h)(n t + h). (3) 

Here n c is the number of particles emitted due to zero and first-order QCD processes 
and n t is the multiplicity of the proton remnant without counting the particles 
migrating from the first-order QCD processes. From (j3|) one obtains 

p=(n c h)-(n c )(h)-(h 2 ) + (h) 2 . (4) 

In this expression, the contributions from the remnant multiplicity n t cancel since 
we consider the case when n t is independent of n c and h. This assumption means 
that the only dominant effect leading to the correlation between n c and n t is the 
first-order QCD migration shown in Fig. 3, rather than non-perturbative effects. 
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The validity of this assumption will be tested using a Monte Carlo simulation in 
Sect. |. 

We do not consider the configuration a) anymore since it cannot contribute to 
the correlation (h = 0). This configuration is rather similar to the QPM events but 
only with large transverse momenta of both jets. Because of this, it is possible that 
some particles migrate to the target region due to the parton showering or resonance 
decays. However, these effects are expected to be small. 

Now let us define the production rate R b for the back-to-back jets (see b) and c) 
in Fig. 3) and the production rate Rf for the events with both jets moving to the 
target region (see d)) as 

where N b is the number of back-to-back jet events, Nf is the number of events with- 
out activity in the current region and N ev is the total number of events. We consider 
these definitions in the limit iV cv — > oo, so that R b and Rf are the probabilities of 
having each configuration. 

For every first-order QCD event with a total multiplicity w, the number h of 
particles moving to the target region is about w/2 for the back-to-back jets, and w 
for events without particles in the current region. Using this estimate, one can write 
the following set of relations: 

(n c h) ~ (w 2 ) [0.5R b + R f ] (6) 



(fi c )(h) ~ (w)(h c ) [O.BR b + R f ] (7) 
(h 2 ) ~ (w 2 ) [0.25/4 + Rf] (8) 

(h) 2 ~ (w) 2 [0.5R b + R f ] 2 . (9) 

Note that the averaging for w is performed only over the relevant first-order QCD 

If 



events. The term (§) can safely be ignored since Rb,Rf << 1, (w) 2 < (w 2 ). Com- 
bining (H)-© together, one obtains from (§) 

p~-(w) (n c ) [0.5R b + R f ] + 0.25 R b (w 2 ). (10) 
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The similarity with e + e~ annihilation makes it possible to estimate (w 2 ). The 
distribution of the multiplicity w in the two-jet events of the e + e~ collisions is usually 
well fitted by a negative-binomial distribution for which the following relation holds 

( w 2 ) = (w) 2 (i + at 1 ) + (w) (li) 

with a free parameter k which is much larger than one (k = oo for a Poisson 
distribution). Using this relation, one obtains 

p~-A b R b -A f Rf (12) 



A b = (w)(n c ) 



1 {w) <l + k-i)- 



(13) 



2 4(n c ) v ' 4(n c ). 

A f = (w)(n c ). (14) 

Below we shall neglect the last term in expression (|13[) which is important only for 
events at small Q 2 when (n c ) < 1. 

The following qualitative predictions can be obtained from ([12]): 

1) For small Q 2 , the multiplicity distribution of the QCD induced events has 
large k, so that it is close to a Poisson distribution. Therefore, for similar values of 
(n c ) and (w), the value of p is negative (anticorrelations). Indeed the correlations 
are negative if 

This relation is expected to be a rather good estimate; From e + e results one expects 
that the QCD induced events well increase the particle multiplicity, however, the 
average first-order QCD multiplicity (w) cannot be much larger than the overall 
average multiplicity (n c ) of all events (including the first-order QCD). 

2) The values of k decreases with increasing Q 2 . This means that p increases 
with energy and even can change sign. The increase with Q 2 is mainly determined 
by the evolution of the average multiplicities (w) and (n c ) as a function of Q 2 . 

3) Let us consider the most interesting case when Q 2 is small and fixed. Since, 
for the current region of the Breit frame, the evolution of multiplicity distribution 
is determined only by Q 2 ||13|| , A b and Af in ( |T2D do not depend on x. Therefore, 
( |T2"D has an x-dependence determined by the production rates R b and Rf. 
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Since the only dominant process at small Q 2 is the BGF, p in ([12]) is mainly 
determined by the behavior of the BGF production rate which increases with de- 
creasing x due to increase of the gluon density inside the proton. Therefore, one can 
expect the magnitude of p increases with decreasing x. 

As we see, one of the most striking features of the Breit frame is a negative value 
of the long-range correlations. Such a prediction is rather unusual for the forward- 
backward correlations studied so far. For the DIS processes in the Breit frame, this 
property is quite clear intuitively: if one or two jets move to the target region, then 
the fewer particles are observed in the current region, the more particles can be 
found in the target region and vice versa. 

Below we shall see that the analytical observations discussed above are in good 
agreement with a Monte Carlo simulation. 



4 Monte Carlo study 

We now illustrate the points discussed above using the LEPTO 6.5 Monte Carlo 



model | 14fl . The model has been tuned as described in |T5|]. The hard process in 
LEPTO is described by a leading order matrix element. Below the matrix-element 
cut-off, parton emission is based on the parton shower described by the Dokshitzer- 



Gribov-Lipatov-Altarelli-Parisi evolution equation. JETSET Monte Carlo |TE[ based 



on the LUND String Fragmentation Model is used to describe hadronization. Al- 



though the production rate of (2+1) jets seems to be underestimated in LEPTO flPfll , 
this Monte Carlo model should be more adequate for illustrating the points discussed 
above since it is based on an exact first-order QCD matrix element calculation. 

To generate DIS events, the energy of the positron and that of the proton is 
chosen to be 27.5 GeV and 820 GeV, respectively. We use the following cuts 

Q 2 > 10 GeV y < 0.95 E > 10 GeV 

where y is the relative energy transfered from the electron to the proton in the 
proton rest frame, E is the energy of scattered electron. In total, 200k events are 
generated for each of the measurements to be discussed below. 
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4.1 First-Order QCD Rates 

Let us first consider the first-order QCD rate for the BGF and QCDC. We define 
the production rates -Rbgf and -Rqcdc for the BGF and QCDC, respectively, as 

R _ ^BGF R _ iVqcDC , 1fi x 

-KBGF - -7T -KQCDC - — l-LOj 

iv ev -'•ev 

where iV B GF and Vqcdc are the numbers of the BGF and QCDC events. Experi- 
mentally, the information about sum of these rates can be obtained from the study 
of (2+1) jet rate -R2+1 [^8j[19[] . Note that -Rbgf + -Rqcdc is not equivalent to -R2+1. 



The latter depends on the jet algorithm and a resolution scale y cut to define jets. 
Note that the normalizations in -R2+1 and flT()| ) can also be different. 

For the Monte Carlo study, ([TJ]) can directly be obtained by counting the first- 
order QCD eventsQ The production rates in LEPTO are derived numerically as the 
integral of the relevant first-order matrix elements. They depend on cut-offs on the 
matrix elements and involve the GRV94 parton density parameterization |2D[ . 



Fig. 4 shows the production rates of BGF and QCDC as a function of Q 2 and x. 
For the latter figure, the cut Q 2 < 50 GeV 2 is used to constrain the effect of increase 
of Q 2 with increasing x. 

For Q 2 variable, the BGF rate rises with increasing Q 2 for Q 2 < 100 GeV 2 , and 
then it falls. For x variable, the BGF rate increases with decreasing x. Since for 
all range of x studied the average value of Q 2 only increases from ~ 19 GeV 2 (for 
(x) ~ 0.0007) to 21 GeV 2 (for (x) ~ 0.02), such a behavior is mainly because of the 
variations of x due to an increase of the gluon density inside proton. 

4.2 Current-Target Multiplicity Correlations 

Fig. 5 shows the behavior of covariance p as a function of Q 2 and x. The values of 
p for the QPM are near zero, i.e. there is no strong linear interdependence between 
the current and target regions^. This illustrates the fact that the LUND String 
Fragmentation Model used by LEPTO does not produce large correlations during 
the formation and an independent breaking of strings stretched between the current- 



region showering partons and the remnants. This is further illustrated in |2T 



1 The parameter LST(24) specifies the type of the first-order QCD event in LEPTO model. 
2 Note again that, generally, a non-linear interdependence may exist. 
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Adding BGF and QCDC events leads to the negative correlations expected. The 
magnitude of p rises with increasing Q 2 (for Q 2 < 100 GeV 2 ) and then falls. Such 
a behavior is due to the similar trend of the BGF rate shown in Fig. 4. Note that 
a decrease of k in ( |i~2"D can also contribute to the decrease of the correlation. For 
Q 2 < 100 GeV 2 , symbols for QPM+BGF and the LEPTO default (QCD) are very 
close to each other. 

The dependence of p as a function of x is mainly due to the variation of x, 
since there are no large variations in Q 2 for the cuts applied. The magnitude of p 
increases with decreasing x. Symbols for the QPM with all first-order QCD effects 
(QPM+BGF+QCDC) and for QPM+BGF events are very close, i.e. the boson- 
gluon fusion is the main source of the correlation. For larger x, symbols become 
distinguishable since the contribution of the BGF becomes smaller at large x. 

According to analytical expression flT2]) , there exists a linear relationship between 
the p and probabilities Rb and Rf determining the BGF production rate. The solid 
line shown in Fig. 5 illustrates the BGF rate (see Fig. 4) multiplied by the scale 
factor —6.5. The shaded band shows the statistical errors for BGF rate. As seen, 
the behavior of p follows that of BGF rate rather well. 

Note that for very small x, one has to expect a deviation from the linear rela- 
tionship between p and BGF rate since the non-linear term omitted in ( |T0"D cannot 
be longer neglected. 

5 Conclusion 

Long-range correlations in DIS were investigated using the Breit frame. It was shown 
that the correlations between the current and target regions of the Breit frame are 
sensitive to the first-order QCD processes leading to the (2+1) jets. In particular, 
at small x and a restricted interval of Q 2 , the strength of the long-range correlations 
is mainly determined by the boson-gluon fusion process. One of the distinguishing 
properties of the correlations is their negative value and rise in the magnitude with 
decreasing x. 

This method can be used to study the production rates of the boson-gluon fusion 
process, to understand better the gluon density inside proton and to discriminate 
between different Monte Carlo models used to simulate the deep-inelastic processes. 
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The approach has the advantage that it is simple and does not involve any jet algo- 
rithm and a resolution scale necessary to determine jet, since it is based on purely 
topological properties of (2+1) jets in the Breit frame. Therefore, the advantage of 
this method is absence of the systematical uncertainties connected with the ambi- 
guity in determining jets. 

Note also that the jet definitions are not able to distinguish between the proton 
remnant and QCD induced jets if they are close to each other. This problem is 
avoided in the proposed jet-rate measurement: The strength of the correlations is 
determined on the basis of an enhancement of the multiplicity in the target region 
for the first-order QCD events, rather than resolving well-separated jets. 

The current-target correlations can be used to study a difference in particle 
spectra between the current region of DIS and a single hemisphere of e + e~ interac- 
tions. In contrast to DIS processes exhibiting anticorrelations in the Breit frame, 
the forward-backward correlations between two opposite rapidity hemispheres of 
e + e~ annihilation are positive This may lead to discrepancies in attempts to 

compare the current region of DIS with a single hemisphere of e + e~ interactions. 

Estimating the current-target multiplicity correlations analytically, we did not 
take into account possible color exchange effects between the fragmentation of the 
outgoing partons and that of the remnant. Also we did not discuss high-order QCD 
corrections. These topics are rather important for a quantitative confrontation of 
the correlations with the data and have to be studied in the future. 
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Figure captions 



• Figure 1: a) Diagram for the neutral current deep- inelastic scattering in the 
QPM; b) A schematic representation of the Breit frame. Particles with p z < 
belong to the current region. Particles with p z > form the target region. 

• Figure 2: Diagrams in the first-order perturbative QCD: a) Boson-Gluon 
Fusion (BGF); b) QCD Compton scattering (QCDC). For the latter diagram 
the gluon can also be radiated before the interaction with 7*. 

• Figure 3: Typical configurations of (2+1) jets in the Breit frame [0. "1" 
denotes the proton remnant, "2" and "3" are jets due to the first-order QCD 
corrections to the parton model. Jets with p z < form the current region. 
For p z > 0, they belong to the target region. 

• Figure 4: Production rates for the BGF and QCDC as a function of Q 2 and 
x. 

• Figure 5: The values of the covariance p for different bins in {Q 2 ) and (x) 
obtained from the LEPTO 6.5 Monte Carlo model. We show four event topolo- 
gies: 1) quark-parton model (QPM); 2) QPM with the first-order QCD events 
(QPM+BGF+QCDC) 3) QPM with BGF events; 4) QPM with QCDC events. 
For the latter figure, symbols for QPM+BGF+QCDC and QPM+BGF are on 
top of each other. 
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Figure 1 
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Figure 2 




Figure 3 
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Figure 5 
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